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micromagnetic technique 



o: 
o. 

(N; 

Oh' 
CD ■ 

on: 
00 : 



C/2 



Vladimir L. Safonov and H. Neal Bertram 

Center for Magnetic Recording Research, University of California 

- San Diego, 

9500 Oilman Dr., La Jolla, CA 92093-0401, U.S.A. 

(February 1, 2008) 



a 



X3 

o 



> 

m 

a^ 
o 
o 
o 



I 

o 
o 



% 



Magnetization reversal in a fine ferromagnetic grain is 
simulated for the case of an instantaneously applied rever- 
sal magnetic field. The Hamiltonian of the system contains 
the exchange interaction, the uniaxial anisotropy, the Zee- 
man energy and the dipole-dipole interactions. A cubic grain 
is discretized into 64 cubic subgrains and the coupled gyro- 
magnetic equations of motion are solved without phenomeno- 
logical damping. A new scheme to solve these equations is 
developed that utilizes only two variables per sub-cube mag- 
netization and strictly conserves the absolute magnitude. The 
initial stage of reversal is uniform rotation followed by a non- 
linear excitation of nonuniform magnetic oscillations driven 
by this uniform mode. An excess of the initial Zeeman energy 
is transformed into nonlinear spin waves, allowing the average 
magnetization to substantially reverse. The process of mag- 
netization reversal in fine quasi-single-domain grain exhibits 
general features of Hamiltonian wave systems with nonlinear 
difi'usion. This nonlinear diffusion is forbidden for either a 
strong reversal field and/or a small grain size. 

75.50.Tt, 75.40.Mg, 76.60.Es 



I. INTRODUCTION 

The problem of magnetization reversal in a fine single- 
domain grain is of great importance in magnetic record- 
ing physics. Characteristic reversal times under strong 
reversal magnetic field give a physical limitation to both 
data rates and system signal to noise ratio. Dynamic 
magnetization reversal depends on the system or mag- 
netic energy relaxation mechanisms. In addition to en- 
ergy absorption to the "thermal bath" of the host lat- 
tice, magnetization reversal can occur by excitation of 
nonlinear 'spin-wave' like modes. In the present paper 
we explore in detail the excitation of these modes in a 
quasi-single-domain particle. This magnetization rever- 
sal process exhibits features of several problems of non- 
linear and solid state physics, such as kinetics of an orien- 
tational phase transition and energy transfer in nonlinear 
systems of oscillators. Therefore the problem is of inter- 
est to fundamental science. 

Dynamic magnetization reversal begins when the ex- 
ternal magnetic field reaches and exceeds the so-called 
'nucleation' field and the original remanent state be- 
comes unstable. Typically these dynamical processes in 
both single and multidomain magnetic systems are solved 



by integration of the Landau-Lifshitz equations M with 
phenomenological damping [^. This form of magneti- 
zation dynamics describes a coherent rotation of each 
discretized subgrain, in which the magnetization magni- 
tude is conserved. Such approximation is fair for small 
excitations, as in linear ferromagnetic resonance P-p[, or 
for some particular model cases when an excess of energy 
directly goes to the thermal bath without excitation of 
spin waves. 

Recently we have explored magnetization reversal in 
quasi-single-domain grain by neglecting phenomenolog- 
ical damping and specifically exploring nonlinear spin- 
wave like magnetic excitations Q, |01- By quasi-single- 
domain we mean grains whose dimension is on the order 
of the domain wall width or in some cases the exchange 
length. Similar analysis has also been done on rever- 
sal in thin films 1^ , H . In all these cases, according to 
Suhl's speculation pof , the switching process should be 
rapid, governed by nonlinear spin-wave excitations. Re- 
cently spin-waves excitated by large magnetization rota- 
tions have been observed experimentally Q . 

By neglecting phenomenological damping, the total 
system magnetic energy is conserved. These simulations 
have demonstrated that the excess of Zeeman energy is 
transformed to the exchange, anisotropy and magneto- 
static energies by nonlinear spin-wave excitations. The 
purpose of this paper is to present a detailed study of 
magnetization reversal simulations focusing on magnetic 
interactions in a single grain. 

In addition, we describe here a new numerical scheme 
of solving the Landau-Lifshitz equations. Typically 
numerical micromagnetic simulations invoke the three 
Cartesian coordinates of each discretization cube. At 
each integration time step the magnetization is modi- 
fied to satisfy the condition of constant magnitude. This 
scheme is preferable to utilization of the natural two- 
parameter polar coordinate scheme because of singulari- 
ties. In polar coordinates there are always two singular- 
ities at the poles. In this new scheme we utilize two 
coordinates per magnetization sub-cube, however here 
these coordinates are similar to linear combinations of 
the transverse magnetization components. Using a sim- 
ple scheme, all singularieties are avoided. 

In section 2 we discuss the micromagnetic model for 
our quasi-single-domain grain. In section 3 we describe 
the implementation of our new numerical model. In sec- 
tions 4, 5 results of simulations are considered in detail. 



Then we discuss obtained results and draw conclusions. 



II. MODEL 

The magnetic energy (Hamiltonian) of the grain can 
be written in the form: 



H — j {Ucxch 
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Here 
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is the exchange energy density. 
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is the uniaxial anisotropy energy density, 
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is the Zeeman energy density and Z^dmag is the magneto- 
static energy density. 

For simplicity we shall consider a micromagnetic grain 
as a system ofn = 4x4x4 (see, Fig.l) coupled cubic 
subgrains of the volume V = L^ , where L is the linear 
size of sub-cube. Each sub-cube can be characterized as 
a classical spin: 



-M,y/ft7, 



(5) 



where M^ is the vector magnetization of j— th sub-cube, 
h is Plank's constant and 7 > is the gyromagnetic ratio. 
|Mj| = Ms, where Mg is the saturation magnetization, 
and therefore, |Sj | = 5 = MsV/fv). Each spin is localized 
in the center of its sub-cube. The characterization of each 
sub-cube magnetization in terms of an effective spin is 
our choice: the problem we consider is purely classical. 
However the spin notation is convenient for the numerical 
scheme presented in Sec. 3. Spin characterization also 
provides a bridge to quantum spin models. 

The Hamiltonian Eq. (nl) of the system in terms of spin 
notation can be written as: 



7i ^ Ti, 
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Here 
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= -^E'(S.-S,v-5^) 



(6) 



(7) 



j¥j' 



describes the exchange interaction between the nearest 
neighbors, J — 2AL/S'^. Note that this form is identical 
to the usual micromagnetic approximation that assumes 



a linear variation of the magnetization between grain cen- 
ters [||. The term —S'^ in (0) yields Hoxch = in the case 
when all spins are oriented in the same direction. 
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is the uniaxial anisotropy energy, where Hk = 2Ku/Ms 
is the anisotropy field. TYanis = corresponds to all spins 
parallel to the anisotropy axis. 



N 



Hz = h-fH • E Sj - h-yHNS 



(9) 



j=i 



is the Zeeman energy. Tiz = if either iJ = or all spins 
are oriented along the external magnetic field direction. 
The magnetostatic interaction of the sub-cubes will be 
approximated by the dipole-dipole interaction energy 
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where Tjj, — Tj — Yj, and Yj , Yj, are the radius- vectors of 

the j and j' spins, respectively. The constant H^J^ is 
chosen so that Tidmag = if all spins are oriented upward 
along the z direction. 

We shall study the problem of time evolution of the 
normalized total magnetization 



1 
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in the system (0) from the initial state when H — and 
the averaged magnetization at i < is oriented "up- 
ward" in the z direction. At the moment t — 0, a strong 
reversal magnetic field H — {H^,0,H^) is applied. The 
component H^ is negative and the transverse component 
i/^ is taken to be relatively small {\H''\ <^ \H^\ ^ Hk)- 
The net vector field is sufficient to give only one energy 
minimum corresponding to an almost "downward" ori- 
entation of the averaged magnetization. Thus the sys- 
tem acquires an excess of Zeeman energy which can be 
later transformed to nonlinear spin excitations contain- 
ing magnetic anisotropy, exchange and dipole-dipole in- 
teraction energies. This transformation will cause the net 
magnetization magnitude |m| to reduce. Mathematically 
the question of energy transfer is described in Appendix 
B. 



III. NUMERICAL SCHEME 

The Landau-Lifshitz equation without damping in 
terms of classical spin (0) can be written as 



TO 



7H, 
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(12) 
(13) 



This equation describe the evolution of coupled spin com- 
ponents S'J, S'J and S^. As time evolves, Eqs.(|l2|),(|l|) 
give two system conservations: the net spin of each sub- 
grain is fixed 



(s^y + {s]y + (s^y = s' 



(14) 



and the total magnetic energy remains constant. Thus 
only two independent variables exist for each j-th site. 

It is natural to write down only two differential equa- 
tions for each spin in terms of spin components or other 
generalized variables. The most popular form of "two- 
variable" Landau-Lifshitz equation is known to be in 
spherical coordinates. However the use of spherical coor- 
dinates for an arbitrary motion of spin gives a stiff system 
of differential equations with singularities at the poles. 
Moreover, numerical calculations of trigonometric func- 
tions typically take much longer time than simple arith- 
metic operations. This is why the spherical coordinate 
approach is not usually used in micromagnetic calcula- 
tions. In practice the direct method of numerical solution 
of th e Lan dau-Lifshitz equation is used (see, for exam- 
ple. 
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Dynamic equations for all three variables 
are integrated with continuum numerical modification so 
that the Eq. ( [l4| ) is maintained. 

In this paper we present a new scheme for numerical 
solution of the Landau-Lifshitz equation ( [12|), in which 
only two variables are utilized and condition (|14[) is auto- 
matically maintained. We introduce the real and dimen- 
sionless variables qj (as a generalized coordinate) and pj 
(as a generalized momentum) as follows 



S! 



qj^{S + SI)/2, 



5f =: 5 - 



(5 + 5|)/2, 
(9?+P?)/2. 



(15) 



These transformations are exact and can be considered 
to be a form of classical Holstein-Primakoff ||l3| transfor- 
mation (see Appendix). Schlomann ||l^ also presented a 
classical spin in a similar form. 

The equations of motion for qj and pj are the exact 
analog of the Landau-Lifshitz equations (]lj) and can be 
written in the form 



(16a) 
(16b) 



dqj 
dT 
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dT 


Hk dq-j 



where a dimcnsionless time r = 7-ffK^ is introduced. Cor- 
responding spin derivatives are 



dsi 

dqj 

dsi 

dpj 

dsi 

dpj 
dSl 

dqj 




(17) 



It is necessary to mention that Eqs.(16a), ( |l6b| ) con- 
tain a singularity at S* -I- 5*1 = 0. To avoid this singu- 



larity we shall use (16a), (16b) only in the hemisphere 
|Sj| = S" with S"! > 0. This is illustrated in Fig.2. Fol- 
lowing Eq.(5), negative magnetization corresponds to a 
positive spin and the {p, q) coordinates are in the upper 
half hemisphere. Initially before reversal, when Mz > 
it is convenient to introduce analogous variables Qj and 
Pj (see Appendix A) for the hemisphere |Sj| = S with 
S-f < : 



S'f^ 



^,? = -P. 
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S+{QUP^)/2 



(18) 



The equations of motion for Qj and Pj have the form: 

(19a) 



dQj 
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dT Hk 

and spin derivatives are 
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(19b) 
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dP, 
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A continuous transition from one set of variable to an- 
other one is defined by 
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Pj 



-P. 
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(21) 



and 



9j 



IS + S^ 



Pj = -Pj 



IS + S^ 



(22) 



Thus, we begin the calculations with Pj{0) and Qj{Q)- 
After a short integration time At we find Pj{At) 
and Qj{At) and calculate corresponding Sj, Sj and 5*1. 
If iS*! < 0, we continue integration of the Landau-Lifshitz 
equations with Pj and Qj. In the case when S^ > 0, 
the j—th spin trajectory crosses the equator and enters 
the upper hemisphere. To avoid the singularity in the 
Pj and Qj coordinates in the upper pole (see Fig. 2), the 
new variables pj and qj are introduced by Eq. (Elf) . Inte- 
gration of the Landau-Lifshitz equations using pj and Qj 
until the z-component of the spin changes sign {S^ < 0). 
At that point we change back to the Pj and Qj coordi- 
nates. 
J^should be noted that the Eqs.(|l6a|), (|l6h|), ( pJa]) and 



( |l9q ) have simple algebraic operations only and therefore 
numerical integration of these equations do not contain 
time consuming operations. 

A characteristic time scale of the system is I/'^Hk- 
The role of exchange interaction is defined by the pa- 
rameter Hc/Hk, where H^ = JS/Hj = 2A/MsL'^ is the 
exchange field from the nearest neighbor. One can write 



i7e/i/K=(n2/Vl6)((5w/(5s)^ 



(23) 



where S^j = 4:{A/Kuy^^ is the domain wall width, 6s = 
ri^/^i = 4L is the linear sample size. The exchange in- 
teraction energy can also be characterized by, so-called, 
"exchange length" (5ex = A^^'^/Mg . Therefore one can 
write 



HJHk = 2n2/3(A#,/i7K)(<5cx/5s)' 



(24) 



We shall focus on the case of relatively strong uniax- 
ial anisotropy when Ms/Hk ^ 1 (which is valid for real 
magnetic recording grains). In the zeroth order of this 
small parameter we can neglect the dipole-dipole interac- 
tion in the system. However as it will be demonstrated in 
the next section that even relatively small dipole-dipole 
interaction, which absorbs very small amount of the Zee- 
man energy plays an important role in the magnetization 
reversal. 



A. Numerical implementation and initial conditions 

The system of Landau-Lifshitz equations has been nu- 
merically solved for n — 64. High-accuracy solutions to 
these ordinary differential equations were obtained with 
minimal computational effort by utilizing the Bulirsch- 
Stoer method with adaptive stepsize control [Q. Inte- 
gration steps were about four orders less than the char- 
acteristic time I/'jHk- 



The equations of motion ( |l6aD , ( |l6h| ), ( |l9aD , ( |l9b| ) 
contain relative parameters only. For simulations these 
parameters were chosen to correspond to real magnetic 
recording media. For example, with A = 10~^ erg/cm, 
Hk = 7 kOe and Mg — 250 Oe, the scaled parameters 
are Sy^ — 43 nm, 6cx = 40 nm and AIs/Hk = 0.036 . 

Nonlinearities can only be excited were there are devi- 
ations in the magnetization state from perfect alignment. 
We include this by allowing the system to initially be in 
thermal equilibrium. We assign randomly a total energy 
deviation of 64^3^. Practically we randomly deviate the 
magnetization of each cell within a fixed limit. This limit 
is adjusted to achieve the above total energy deviation. 
To obtain a true thermal equilibrium the Landau-Lifshitz 
equations are solved, given the above initialization, for a 
time ^HyJ^q ~ 200. At this point the system is prepared 
for the application of the reversal field. 



IV. MAGNETIZATION REVERSAL WITHOUT 
DIPOLE INTERACTIONS 

First we shall consider the case when the dipole-dipole 
interactions in 7Hcftj ( [l3| ) are neglected. Fig.3a,b show 
the time dependencies of mz{t) and |m(i)| in the case 
when the reverse magnetic field has a very small devia- 
tion from the anisotropy axis {Hx/H^ = 0.02). Fig. 3b is 
simply an extension of the time axis of Fig. 3a. The field 
magnitude is slightly less than the nucleation field for uni- 
form rotation {Hnuc — — O.OiJx)- This field corresponds 
(by numerical test) to the lowest reverse field required 
to produce nonlinear excitations. The sample dimen- 
sion is 16% and 25% greater than the domain wall width 
and exchange length, respectively. Initially the longitu- 
dinal component of magnetization decreases monotoni- 
cally from niz ~ 1 to ~ 0.7 and then at a 'critical' time 
^H}<itc — 70, oscillations of rriz occur. In contrast, the 
net magnitude of magnetization remains practically con- 
stant |m| ~ 1 until tc and then abruptly decreases ex- 
hibiting irregular oscillations. For times less than tc the 
average magnetization simply precesses about the reverse 
field. For times greater tc the excitation of nonlinear spin 
waves becomes pronounced. These excitations both ab- 
sorb Zeeman energy and reduce the magnetization. 

This behavior is typical of nonlinear spin-wave pro- 
cesses whose amplitudes grow exponentially from an ini- 
tial thermal level. The decrease in net magnetization 
becomes observable at virtually a critical time tc when 
the magnetic moment of the (transverse) excited waves 
becomes comparable with the net magnetic moment. At 
the same time the system nonlinearities restrict the ex- 
ponential growth. The oscillations in rriz that occur for 
t > tc primarily correspond to the uniform rotation: the 
frequency increases with decreasing rriz . If spin- wave ex- 
citations did not occur, oscillations in tUx, my, rriz would 
reflect uniform rotation only and at every instant sum to 
give constant |in| ~ 1. 



With increasing time the average niz continually de- 
creases reaching, for this example an asymptotic level 
in the vicinity of m^ = (at least for times less than 
"fHyit ^ 3000). There is no reversal of mz{t). The aver- 
age magnetization |ni(i)| also decreases by about « 60%. 
As time proceeds, the spin-wave oscillations exhibit in- 
creased chaotic behavior. 

Fig. 4 shows similar time dependences as in Fig. 3 but 
for the case of increased angle of the reverse magnetic 
field from the anisotropy axis {Hx/ H^ — 0.10). The field 
magnitude {Hz/Hk — —0.70) also corresponds to the 
onset of spin-wave excitation. As expected, the magne- 
tization dynamics is much faster in this case. We can 
see the beginning of nonlinear spin-wave excitation at 
7-ffK^c — 30, where the ruz has reduced due to pure rota- 
tion to ~ 0.3. The decrease of longitudinal magnetization 
is more substantial, passing through zero at jHjito ~ 75 
and reaching about w 55% of its nominal value in the 
opposite direction after 7-ffKi ^ 1000. Similar to the 
case shown in Fig. 3, the average magnetization |m(i)| 
decreases to « 60%. Note that the amplitudes of chaotic 
oscillations of |ni(i)| are relatively smaller than in the 
case of Fig. 3 (compare the long time behavior of Fig. 3b 
and Fig. 4b). Here most of the Zeeman energy is trans- 
formed into high frequency exchange-dominant spin- wave 
modes, which have relatively small magnetic moments. 

As discussed earlier, the total magnetic energy remains 
constant in these dynamic processes. In Fig. 5a, b we plot 
the evolution of the Zeeman, exchange and anisotropy 
energies (Fig. 5b is an extended scale of Fig. 5a). The re- 
verse angle is intermediate between Fig. 3 and 4. Initially, 
for times less than tc (in this case ^HKtc — 45), uniform 
rotation yields a decrease of the Zeeman energy and an 
increase of the anisotropy energy. At the onset of nonlin- 
ear spin-wave excitation we see, in addition, an increase 
of the exchange energy. The wave-like nature of these 
excitations can be seen in Fig. 5a in the beating of the 
exchange and anisotropy energies. The nonlinear excita- 
tion of chaotic oscillations of the magnetization can be 
thought of as an "overheating of the magnetic system" . 
The more exchange magnetic modes that are involved 
in the process, the smaller the averaged magnetization 
|m(i)|. 

Fig.6a-d demonstrate the time dependencies of mz{t) 
and |m(t)| for various magnetic fields: — 0.78iJK, 
-0.79iJK, -0.90ifK and -I.OHk, respectively. As in 
Fig. 5 the field angle is fixed at Hx/Hz — 0.05 and rel- 
ative sizes are unchanged. For this angle the Stoncr- 
Wohlfarth nucleation field is — 0.83iJK- In this case 
the onset field for nonlinear spin-wave excitation is Hz 
=—0.79Hk- There is no reversal below this field as shown 
in Fig. 6a. The magnetization evolution at the onset field 
is shown in Fig. 6b (similar to Figs. 3, 4). Magnetization 
dynamics for larger fields are shown in Figs. 6c, d. A no- 
table features can be seen: the critical time for nonlinear 
spin-wave excitation increases with increasing field. The 
increase of tc allows for several cycles of uniform rotation 
to occur prior to the nonlinear spin-wave excitation. No 



reversal dynamics was observed beyond Hz < —I.IHk 
(in this particular case). Essentially as the reverse field 
is increased beyond the spin-wave onset field, tc rapidly 
increases and becomes infinite. These phenomena also 
occur for other field angles {Hx/Hz). 

As indicated in Figs.6a-d the asymptotic, long time 
value of niz increases with increasing applied field (i.e. 
less magnetization reversal occurs). In Fig. 7 we explic- 
itly plot the asymptotic values of averaged {mz{t)) and 
(|m(i)|) versus apphed field. It should be noted that the 
excitation of nonlinear spin waves occurs only for a nar- 
row band of applied fields. 

Fig. 8 shows the evolution of magnetization at various 
sample sizes. We see that a decrease of the sample size 
relative to the domain wall width leads to a suppression 
of magnetization reversal. This fact is also can be in- 
terpreted in terms of nonlinear spin- wave processes. By 
decreasing the sample size, we increase the frequencies of 
the non-uniform spin-wave modes and therefore restrict 
the number of resonance nonlinear processes that can 
contribute to a decay of the uniform mode. Analysis for 
larger sample would require finer discretization. 



V. MAGNETIZATION REVERSAL WITH 
DIPOLE INTERACTIONS 

In this section we show numerical simulations that in- 
clude the dipole-dipole interactions in 7Hcffj ([li^). The 
relative role of dipole-dipole interaction is small because 
we assume a small magnetostatic parameter Ms/ Hk <C 
1. In this case we can not expect a noticeable absorption 
of the Zeeman energy solely by the dipole-dipole reser- 
voir. Nevertheless, the dipole-dipole interaction plays an 
important role in the process of energy transfer: it re- 
duces the symmetry of the system and therefore opens 
new channels of nonlinear spin-wave interactions. 

Let us compare the results of Fig. 3 obtained without 
dipole interaction and corresponding results in Fig. 9 with 
the dipole interaction. The onset time tc of nonlinear 
excitation is about the same. However mz{t) decreases 
more rapidly, reverses and reaches an asymptotic value of 
~ —0.5. This is in contrast to Fig. 2 where the asymptotic 
mz{t) does not reverse. In this case the average magneti- 
zation |m(t)| decreases to ~ 50%, slightly more than the 
w 60% in Fig. 3. The amplitudes of chaotic oscillations of 
|m(t)| in this case are much smaller than occur in Fig. 3. 
The excess Zeeman energy is transformed into exchange 
modes, which have smaller magnetic moments. 

Plotted in Fig. 10 is the magnetization dynamics for a 
larger field angle {H^/Hz = 0.10). As compared with 
Fig. 4, including magnetostatic interactions slightly in- 
creases the reversal rate. The longitudinal magnetization 
passes through zero at ^HyJ^q ~ 55 and, similar to Fig. 4, 
reaches about w 55% of its nominal value in the opposite 
direction after jHKt > 500. The average magnetization 
|m(i)| is decreased to « 60%. 



Fig. 5c, d demonstrate typical transformations of Zee- 
man, anisotropy and exchange energies for short (c) and 
long (d) time intervals for magnetization dynamics with 
dipole interactions. We see faster dynamics of energy 
transformations than shown in Fig. 5a, b. It should be 
noted that the energy absorbed by dipole interactions is 
relatively small (within ±0.01 in relative units) as ex- 
pected and therefore not shown. However the magneto- 
static interactions change the balance the exchange and 
anisotropy energies. 

In Fig. 11 are shown the time dependencies mz{t) and 
|m(i) I for various reversal fields. In agreement with Fig. 6, 
the reversal dynamics is suppressed with increasing the 
strength of reversal field but the range of reversal dy- 
namics is much broader. We did not observe any rever- 
sal dynamics at H^ < — 3-ffK in this particular case. We 
believe that the 4x4x4 discretization is not a major 
limiting factor. For the case of n = 1000 subgrains no 
reversal dynamics was observed at H^ < —AHy^ 0. 

In Fig. 12 we plot the asymptotic relative change in 
Zeeman energy Ai?z = (1 — (mz))(— if^/i?K) versus 
the applied magnetic field for the parameters of Fig. 11 
(solid dots). Even though there is less asymptotic mag- 
netization decrease as the field is increased, the rela- 
tive change in Zeeman energy does increase. After the 
onset of nonlinear spin-wave excitation, in this case at 
Hz — — 0.79iJK, the amount of spin-wave energy that 
can be excited increases to a maximum value. In addi- 
tion (open dots) , the change in Zeeman energy is plotted 
corresponding to Fig. 5. for the case of no magnetostatic 
interaction. The same trend occurs but only over a small 
interval of field values. 

Fig. 13 shows the evolution of magnetization at various 
sample sizes. As in Fig. 8, the decrease of the sample size 
relative to the domain wall width leads to a suppression 
of magnetization reversal. The minimum size for spin- 
wave excitation remains the same, however differences 
are seen in the amplitude of the chaotic excitations. 

We have also examined the effect of grain size relative 
to the exchange length at fixed relative domain wall size 
((5s/(5w — 1.16). The system evolution is more sensitive 
to the change of exchange length if the dipole-dipole in- 
teraction is taken into account. The primary effect is 
that with a smaller relative exchange length the rate of 
decrease of the average magnetization is faster. As dis- 
cussed above, the presence of magnetostatic interactions 
allows the excitation of additional spin- wave modes that 
were otherwise degenerate. These modes are more eas- 
ily excited if the sample size is large compared to the 
exchange length. 



VI. DISCUSSION 

Magnetization reversal occurs by reducing the Zeeman 
energy as the magnetization goes from an initial high en- 
ergy state to one of low energy (approximately in the field 



direction). Typically this dynamic process is analyzed 
using the Landau-Lifshitz equation with an appropriate 
phenomenological damping term. Here we have excluded 
phenomenological damping and shown that substantial, 
rapid reversal can occur. In this case the total mag- 
netic energy is conserved, however subject to constraints 
of sample size and field range, the Zeeman energy is re- 
duced by the excitation of nonlinear spin-waves. 

In the present study we have specifically focused on 
a cubic grain sufficiently small to be almost single do- 
main. In this case after uniform rotation over a short 
time period, nonlinear spin waves are excited and the 
average grain magnetization decreases. For the sample 
sizes examined here the reduction was from positive sat- 
uration {-\-Ms) to at most, — 0.6Ms. Due to numerical 
limitations we did not examine large grains and it is pos- 
sible that, with sufficient discretization, almost complete 
reversal might occur for samples much larger than the 
domain wall width. In this study we considered only 
samples with small relative magnetostatic energy con- 
tribution. Studies have been performed on high moment 
low anisotropy thin films, where a large out of plane mag- 
netization reversal torque occurs f^], @. In these cases 
virtually complete reversal almost always occurs with no 
phenomenological damping. 

The decrease and suppression of magnetization rever- 
sal with increasing the reversal field strength initially 
seem unusual. In the discussion of Fig. 11 a limit is sug- 
gested as due to the saturation of the spin-wave energy. 
However, the problem we consider is a general system 
of nonlinear oscillators (nonlinear spin waves). In the 
oscillator problem one mode (here corresponding to co- 
herent rotation of the magnetization) is initially strongly 
excited. One might expect that the system would move 
towards its most probable state (thermal equilibrium) by 
a chaotic mode mixing. With any large reverse field we 
would expect reversal to a new equilibrium state. How- 
ever, the observed reversal decreases and suppression can 
be qualitatively understood as a result of nonlinear spin- 
wave interactions. 

The processes of nonlinear spin- wave decay are defined 
by several factors. First, it is necessary that the fre- 
quency of the initial mode and some portion of the non- 
linear spin-wave spectrum overlap. Second, we require 
sufficiently strong nonlinear interwave coupling. These 
factors determine a threshold amplitude of the decaying 
wave (as an effective pumping field). Third, the ampli- 
tude of initial (decaying) wave must exceed a threshold 
value. When the reversal field strength is increased, all 
these three requirements eventually are not met. As sug- 
gested in Ref. jl^l, increasing the reverse field shifts the 
spin-wave spectrum away from the uniform precession 
mode. Also as the reversal field is increased nonlinear 
coupling decreases. At a certain reversal field the thresh- 
old of nonlinear decay will exceed the amplitude of coher- 
ent magnetization rotation and the reversal process will 
be suppressed. Generally speaking, such a suppression of 
the energy transfer from the excited mode to other modes 



is a common feature of the theory of Hamiltonian systems 
of nonhnear osciUators (see, e.g., ffT))- The most simple 
and popular example of nonlinear coupled oscillators is 
the Fermi-Pasta- Ulam chain (see, e.g., |18f). 

The numerical scheme we have presented is valuable 
here primarily because the number of coordinates is re- 
duced by 2/3. Thus, the computational storage require- 
ment is correspondingly reduced. However, there is an 
additional advantage to this formulation. These coordi- 
nates correspond naturally to nonlinear oscillators that 
describe spin-wave motions. Recently, the problem of 
thermally agitated reversal has been analyzed in a mi- 
cromagnetic grain. Using this nonlinear oscillator repre- 
sentation a simple exact solution has been obtained for 
the first passage time 111, H, (21). 



VII. CONCLUSION 

Magnetization reversal in a quasi-single-domain grain 
can occur by nonlinear spin-wave excitation. In this pro- 
cess an excess of Zeeman energy is transfered to uni- 
axial anisotropy and exchange energies. The dipole- 
dipole interaction plays an important role even if small 
{Ms/Hk <C 1). The process of magnetization reversal in 
a fine grain exhibits general features of system of inter- 
acting oscillators: energy transfer by nonlinear resonance 
processes. This nonlinear process (diffusion) is forbidden 
for either a strong reversal field and/or a small grain size. 

We have presented a new scheme to perform numer- 
ical integration of the Landau-Lifshitz equations. This 
new scheme utilizes only two variables per discretiza- 
tion cell, but strictly conserves the magnetization mag- 
nitude. Comparison with conventional techniques show 
that both schemes require approximately the same com- 
putation time at the same level of accuracy (the new 
scheme takes a slightly shorter time). The main advan- 
tage of the new scheme is that it utilizes a smaller number 
of variables (2/3 of the total magnetization components) 
and therefore requires less memory. 
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APPENDIX A: EQUATIONS OF MOTIONS AND 
CANONICAL VARIABLES 

1. Classical commutations 

Hamilton's equations can be written in a general form 

(Al) 



^n'§^lo,n] 



Here is a function of interest and | ....... ] = 

i{. . . , . . .} denotes a classical analog of commutator and 
{...,...} are the Poisson brackets. Plank's constant h is 



used as a dimensional constant to simplify a dimensional 
structure of the equation and for simple correspondence 
with quantum mechanics. 

Taking into account the Mermin's [Q formula for the 
Poisson brackets for a system of classical spins S^, one 
can write: 



M.«i^.i:s,^^«S 



dSj dSj 



(A2) 



It is simple to show that the commutation rules for the 
classical spins are the same as in quantum mechanics: 

lS],Sf^]^±Sfd,,„ IS+,SJ-,1^2S]S,„. (A3) 
where 5 = S*^ ± iS^ are the circular spin components. 



Substituting O — Sj to Eq. (Al) and taking into ac- 
count that 

.as ^ dn 
= -^as-^"as 

one can obtain the Landau-Lifshitz equation (|l2|). 



2. Holstein-Primakoff transformations 



It may be easily shown that the commutations (A3) 
are valid if 



S+ = a^/sTS^, 



S- = a* VS + S'. 



(A4) 



These relations describe the classical spin realization in 
terms of complex variables a* and a corresponding to 
quantum creation and annihilation operators, respec- 
tively. The transformation ( A4 ) in the quantum case was 
introduced by Holstein and Primakoff [H . Suhl |23] first 
utilized approximate Holstein-Primakoff transformation 
with complex variables to describe magnetic resonance. 
The classical analog of the commutator for a* and a 
is 



dA dB 
da da* 



dBdA 
da da* 



lABl 

The Eq.( |Al|) with (A5) can be rewritten as 
dn/h dS 



da 



as da' 



= -ijUcB ■ 



dS_ 
d^ 



(A5) 



(A6) 



Another form of Holstein-Primakoff transformation 
can be written as 



S+ = b* VS^^5^, 



S- =bVS^^S^. (A7) 

The classical analog of commutator for b* and b is 

dAd^_dBdA 
^ ' ^ - db db* db db* 



and the corresponding equation of motion has the form 

db dn/h as ^^ as 



(A8) 



It should be noted that Eqs.(A6) and ( |Aq ) represent 
other forms of the Landau-Lifshitz equation (n2h . 



3. Real canonical variables 



Let us introduce real variables q and p as 



q + ip 



V2 



Thus the Eqs.(A4) become 



S- = q^{S + S^)/2, 

sy - p ^{s + s-)/2, 

S'- ^S-{q^+p')/2. 



(A9) 



These formulas are exact and can be considered as a form 
of classical Holstein-Primakoff transformation. 
The classical commutator for q and p is: 



^ dq dp dq dp 



(AlO) 



The equations of motion for real and dimensionless q and 
p are 

dq _ dn/h _ dn/h as 

dt dp as dp ' 

dp _ dn/h _ dn/h as 
di ^" dq ~ as~' a^' 



(All) 



We shall use the transformation (A9) in the case, when 



5| > 0. 



In order to exclude stiff equations, for S"! < it is con- 
venient to introduce other variables Q and P as follows 



b = 



Q + iP 



b* = 



Q-iP 



V2 ' V2 

The equations of motion for Q and P have the form: 



dQ _ dn/h as _ as 

d7 " ^s~ ■ dP ~ ^ "'^ ■ ap' 
dP _ dn/h as__ as 
d7~ as~'ag~"^ ''^'dQ' 



(A12) 



APPENDIX B: ENERGY TRANSFER 

Mathematically the question of energy transfer be- 
tween Tiz, "^cxch, Wanis and Hdd- IS associated with the 



Eq. (Al) with = Hz: 



ih^ = inz,nz 



'Lcxch ~r 'Lai 



^dd] 



(Bl) 



It is obvious that |7iz,Wz] — 0. Direct energy transfer 
from the Zeeman energy to the exchange energy is always 
forbidden, because 



[Hz,^: 



cxch 



= 0. 



(B2) 



The commutator | Tiz , Tianis ] depends on the orienta- 
tion of the external magnetic field. If the external mag- 
netic field is parallel to z axis (H = {0,0, H)) then this 
commutator is equal to zero. However, even a slight devi- 
ation of H from the z axis permits a possibility to trans- 
form in sequence the Zeeman energy to the anisotropy 
energy and the latter to the exchange energy as long as 



['^cxchj'^anisl 7^ 0. 



(B3) 



Captions 

Fig.l. 

Model of micromagnetic grain (4x4x4 subgrains). 

Fig.2. 

New coordinates. 

Fig.3. 

Short (a) and long (b) time evolution of magnetization 
in the case of very small deviation of reversal magnetic 
field from the anisotropy axis. 

Fig.4. 

Short (a) and long (b) time evolution of magnetization 
in the case of moderate deviation of reversal magnetic 
field from the anisotropy axis. 

Fig.5. 

Transformations of Zeeman, anisotropy and exchange 
energies for short (a,c) and long (b,d) time intervals. The 
dipole interactions are neglected in (a,b) and are taken 
into account in (c,d). Hz = —0.79Hk is the onset of 
reversal, H^/H^ ^ 0.05, S^/Se^ = 1.08. 

Fig.6. 

Time evolution of magnetization for various reversal 
fields at Hx/Hz — 0.05, 6v,/6ex = 1-08 (dipole interac- 
tions are neglected). 

Fig.7. 

Average asymptotic longitudinal component and abso- 
lute value of magnetization versus applied reverse field. 
Conditions correspond to Fig.6. 

Fig.8 

Time evolution of magnetization for various grain sizes 
Ss/6^. Hz/Hk - -0.79, H,/Hz = 0.05, Jw/^ex = 1.08 
(dipole interactions are neglected). 

Fig.9. 

Short (a) and long (b) time evolution of magnetization 
in the case of very small deviation of reversal magnetic 
field from the anisotropy axis. The strength of reversal 
field corresponds to the onset of reversal (dipole interac- 
tions are taken into account). 

Fig.lO. 

Short (a) and long (b) time evolution of magnetization 
in the case of moderate deviation of reversal magnetic 
field from the anisotropy axis. The strength of reversal 
field corresponds to the onset of reversal (dipole interac- 
tions are taken into account). 

Fig.ll. 

Time evolution of magnetization for various reversal 
fields at Hx/Hz — 0.05, ^w/<5ox = 1-08 (dipole interac- 
tions are taken into account). 

Fig.l2. 

Relative change in Zeeman energy versus applied re- 
versal magnetic field. 

Fig. 13. Time evolution of magnetization for various 
grain sizes 5J5^. Hz/Hk ^ -0.79, H^/Hz = 0.05, 
^w/^ex — 1-08 (dipole interactions are taken into ac- 
count). 
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